On Anisotropy of Ultra-High Energy Cosmic-Rays 

Tamar Kashti 

Weizmann Institute of Science, Israel 



Abstract 

We briefly summarize our study on anisotropy of Ultra-High Energy Cosmic- Rays (UHECRs), 
in which we define a statistics that measures the correlation between UHECRs and Large Scale 
Structure (LSS). We also comment here on recently published paper by Koers and Tinyakov 
that compared our statistics to improved KS statistics. 

Key words: ultra high energy cosmic rays, cosmic rays, large scale structure 
PACS: 95.85.Ry, 98.65.-r 



The origin of cosmic rays of energies > 10 19 eV is a puzzle [11213] . The arrival di- 
rection of UHECRs show no correlation with the galactic disk, which point towards 
extra-galactic origan. A suppression of the spectrum at ~ 5 x 10 19 eV is expected due 
to the GZK suppression |4|5j . The suppression was observed by HIRES [6] and the new 
Auger Observatory [7], see fig.[TJa). Thus, cosmic-rays with energies above ~5x 10 19 eV 
can reach us only from sources with distance below ~ 100 Mpc. In these distances, the 
Universe is not isotropic, implying that correlation with LSS can give us information on 
the origin of UHECRs. 

There are few leading candidates as the sources of UHECRs. Assuming that UHECRs 
are protons, the large magnetic fields needed for acceleration requires source luminosity 
of L > 10 12 L Q r 2 //3 for protons with 10 20 eV, where L G is the Sun luminosity and T 
is Lorentz factor of the magnetized plasma [3]. Therefore, the only known astrophysical 
sources that reaches these energies are Active Galactic Nuclei (AGNs) and sources of 
Gamma-ray bursts. The anisotropy of UHECRs from these sources should be correlated 
with LSS. Another possible source of UHECRs is the decay of new heavy particles coming 
from top-down models pQ, which predict isotropic signal. 

In [8], we derive the expected all sky angular distribution of the UHECR intensity, 
and defined a statistics Xc,ub that measures the correlation between the predicted and 
observed UHECR arrival direction distribution. Following [9] , we consider a model where 
the UHECR flux is produced by cosmological sources of protons tracing the large scale 
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Fig. 1. (a) The spectrum of HIRES and Auger after a shift of AE/E = +23% in the calibration of the 
absolute energy scale of the Auger experiment. The solid line is the spectrum that would be generated 
by a cosmological distribution of sources of protons, with intrinsic spectrum d log n/d log E = —2. (b) 
The positions of the 27 Auger events with energy exceeding 5.7 X 10 19 eV, overlaid on the intensity map 
obtained in the biased model, in galactic coordinates. 

galaxy distribution. We assume that the sources are intrinsically identical and that the 
number density of sources is drawn from a Poisson distribution with an average given by 
6[(S]s(z), where s(z) is the average comoving number density of sources at redshift z and 
b is some bias functional of the local fractional galaxy over density, 8 = 8p/p. The LSS 
galaxy density field is derived from the PSCz catalogue [10]. For the bias functional we 
consider three models: an isotropic (I) model, b[8] = 1; an unbiased (UB) model where 
the source distribution traces the galaxy distribution with b[6] = 1 + 8; and a biased (B) 
model, b[8] = 1 + 8 for 8 > and b[8] = otherwise. 

The statistics we defined in j8] , which measures the correlation between predicted and 
observed UHECR arrival direction distributions, is: 



Here {i} is a set of angular bins, iVj is the number of events detected in bin i, and A^m is 
the average number of events expected to be detected in the M (e.g. isotropic, unbiased, 
biased) model. In order to avoid sensitivity to the possible distortion of the UHECR inten- 
sity map by magnetic fields, we used 6° x 6° angular bins and excluded the Galactic plane 
region, |6| < 12°. The value of Xcub can be straight forwardly calculated using the nu- 
merical representations of the UHECR maps at http : //www. weizmann. ac . il/~waxman/ criso 

The Xq statistics is more sensitive to expected anisotropy signature than the com- 
monly used power spectrum, Ci = Em=— « a im ( e -§- ED; an d the two-point corre- 
lation function, W(D) = ^2j<i ®{P ~ ^ij) ( e -g- [12] )■ This can be seen from table 
[lja). The anisotropy signal is stronger at lower energy: Although the contrast of the 
fluctuations in the UHECR intensity is higher at high energy, the signal becomes weaker 
at higher energies since the number of observed UHECR drops rapidly with energy, see 
table QJb) . In [8] we also show that a few fold increase of the Auger exposure is likely 
to increase the significance to > 99% CL, but not to > 99.9% CL (unless the UHECR 
source density is comparable or larger than that of galaxies) . 
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E > 40 EeV, 100 events 

Xc,ub 
X W ({D = 10,20,30,40}) 

Xc({e = 2}) 



P(I/UB) P(I/B) P(UB/B) 
23% 79% 42% 
7% 12% 10% 

6% 8% 7% 



Xc.UB 
E > 20 EeV (1205 events) 
E > 40 EeV 300 events 
E > 60 EeV (94 events) 
E > 80 EeV (31 events) 



P(I/UB) P(I/B) P(UB/B) 



45% 


94% 


45% 


39% 


94% 


52% 


31% 


87% 


42% 


22% 


63% 


24% 



Table 1 

Probabilities P(M\\M2) to rule out model M\ at a certain confidence level (CL), assuming that the 
UHECR source distribution follows model Mq.. Calculated from 10000 Monte Carlo realizations for Auger 
exposure, (a) Comparing different statistics, Xc,UBi Xc{£) = slf (Qf- ~ Ciso,£) 2 where <r| is the 
variance, and Xw (defined similarly for the two-point correlation function), (b) Comparing different 
energy thresholds with Xc ub statistics. 

Auger reported in [13] a correlation between the arrival direction distribution of 27 
UHECRs of > 5.7 x 10 19 eV and between the angular distribution of low-luminosity 
AGNs included in the V-C AGN catalog [T3J. However, the V-C AGN catalog is merely a 
compilation of AGN data available in the literature, and is therefore incomplete both in 
its sky coverage and in its luminosity coverage, hence the results of [13j are unclear. Using 
the Xc,ub statistics, we analyzed [13] data, see fig. [ljb). According to our analysis, the 
data is inconsistent with isotropy at ~ 98% CL, and consistent with a source distribution 
that traces galaxy density. Note, however, that the optimization of the energy threshold 
made in [13] raises the concern that the significance with which isotropy is ruled out 
maybe overestimated. 

Recently, [15] compared the Xc,ub statistics to an improved KS statistics. They 
showed that for most energy thresholds the Xc,ub statistics have higher probability 
to rule out isotropy (see their table I). However, they get probabilities that are different 
then ours for Xc.ub- This is probably due to the fact that they use a different galaxy 
redshift catalog, the catalog of Kalashev et al. [16], deduced from 2MASS XSC [17]. We 
believe that this catalog is not suitable for the analysis of CR anisotropy, since it includes 
mainly photometric redshifts, which have systematic uncertainty of 30% that depends on 
the luminosity, thus distorting the density field. Note that the future "2MASS Redshift 
Survey" [T7] will improve PSCz, by measuring redshifts of 100,000 galaxies till 85 Mpc. 

References 

[I] P. Bhattacharjee and G. Sigl, Phys. Rept. 327, 109 (2000). 

[2] M. Nagano and A. A. Watson, Rev. Mod. Phys. 72, 689 (2000). 

[3] E. Waxman, Pramana 62, 483 (2004), and his talk in these proceedings. 

[4] K. Greisen, Phys. Rev. Lett. 16, 748 (1966). 

[5] G. T. Zatsepin and V. A. Kuzmin, JETP Lett. 4, 78 (1966). 

[6] R. Abbasi et al [HiRes Collaboration], Phys. Rev. Lett. 100, 101101 (2008). 

[7] J. Abraham et al. [Pierre Auger Collaboration], Phys. Rev. Lett. 101, 061101 (2008). 

[8] T. Kashti and E. Waxman, JCAP 0805, 006 (2008). 

[9] E. Waxman, K. B. Fisher and T. Piran, Astrophys. J. 483, 1 (1997). 

[10] W. Saunders et al, Mon. Not. Roy. Astron. Soc. 317, 55 (2000). 

[II] G. Sigl, F. Miniati and T. A. Ensslin, Phys. Rev. D 70, 043007 (2004). 
[12] D. De Marco, P. Blasi, A. V. Olinto, JCAP 01, 002 (2006). 

[13] Auger Collaboration, Science, 318 (5852), 938; J. Abraham et a/.,Astropart. Phys. 29, 188 (2008). 

[14] M.-P.Veron-Cetty & P.Veron, A&A 455, 773 (2006). 

[15] H. B. J. Koers and P. Tinyakov, arXiv:0812.0860 [astro-ph]. 

[16] O. E. Kalashev, B. A. Khrenov, P. Klimov, S. Sharakin and S. V. Troitsky, J CAP 0803, 003 (2008). 
[17] Jarrett, T.H. 2004, PASA, 21, 396, |http: //www, cf a .harvard . edu/~huchra/2mass/| 



3 



